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Abstract: Interference is a serious cause of performance degradation for IEEE802.15.4 
devices. The effect of concurrent transmissions in IEEE 802.15.4 has been generally 
investigated by means of simulation or experimental activities. In this paper, a mathematical 
framework for the derivation of chip, symbol and packet error probability of a typical 
IEEE 802.15.4 receiver in the presence of interference is proposed. Both non-coherent and 
coherent demodulation schemes are considered by our model under the assumption of the 
absence of thermal noise. Simulation results are also added to assess the validity of the 
mathematical framework when the effect of thermal noise cannot be neglected. Numerical 
results show that the proposed analysis is in agreement with the measurement results on the 
literature under realistic working conditions. 

Keywords: concurrent transmissions; capture effect; IEEE 802.15.4 



1. Introduction 

Among the possible medium access control (MAC) techniques for wireless communication systems, 
the simplicity of random access schemes {i.e., ALOHA, carrier sense multiple access (CSMA)) make 
them suited to be implemented in several standards for short range apphcations [1,2]. Even though 



Sensors 2014, 14 



5623 



mitigation methods can be introduced in random access MAC {i.e., carrier sense multiple access with 
collision avoidance (CSMA-CA)), collisions are not completely avoidable. Nevertheless, some receivers 
have the ability to correctly receive a signal despite a significant level of co-channel interference, and 
collisions do not always lead to a total loss of the collided packets. This co-channel interference tolerance 
is called capture effect [3]. In the presence of concurrent transmissions at the same carrier frequency 
(collisions), packet capture may happen even for low values of the signal-to-interference ratio (SIR). 

The first papers in the literature about the capture effect mostly consider Frequency Modulation (FM) 
demodulators [3-5]. Later, the capture effect has been also studied in a variety of transceivers and 
MAC schemes, including ALOHA networks [6-8], IEEE 802.11 devices [9-11], Bluetooth radios [12] 
and cellular systems [13]. Focusing on IEEE802.15.4, several papers describing experimental results 
can be found in the literature [14-22]. To give some examples, in [14], the packet capture probability 
of a Chipcom CCIOOO transceiver [15] is measured with the aim of exploiting the capture effect for 
collision detection and recovery. Another study [16], carried out again with CCIOOO transceivers, 
which work in the sub-1 GHz band, obtained a threshold for the capture probability for the case of 
one interferer, but unstable results for the multiple interferers. This early work seems to suggest that 
the number of interferers might have an important effect on the capture probability with CCIOOO 
devices. However, in contrast to the previous CCIOOO measurements, successive studies, carried 
out with Chipcon CC2420 transceivers [17], which operate in the 2.4 GHz industrial, scientific and 
medical (ISM) band, show that the capture effect depends on the total interfering power, but it is 
independent of the number of interferers [18,21]; therefore, in our mathematical analysis, we assume 
that the performance is independent of the number of interferers for the 2.4 GHz physical layer (PHY) 
of IEEE 802. 15.4, and as a result, we use only one interfering signal for the mathematical evaluation. A 
behavior on the packet capture similar to that reported in [16,18] is also observed in [20] with Freescale 
MC1224 transceivers [22], which, again, operate in 2.4 GHz. The experiments conducted in interferer 
power-dominant (with respect to the noise) environments in [16,18-20] show a couple of common 
behaviors: (i) the receiver starts capturing the useful packet when signal-to-interference-plus-noise 
ratio (SINR) goes beyond 0 dB; (ii) packet reception rate (PRR) reaches one for values of SINR larger 
than 4 dB (please see in [16], Figure 5 and Figure 16, for CCIOOO and CC2420, respectively; in [18], 
Figure 4 for CC2420; in [19], Figure 7c for CC2420; and in [20], Figure 3 for MC13192). 

Although the experimental results in [16,18-20] agree with each other, there is no theoretical model, 
purely based on mathematical analysis, which can be applied to the 2.4 GHz PHY of IEEE 802.15.4 
to explain such characteristics. Motivated by this consideration, we propose an analytical framework 
to investigate the behavior of the IEEE 802.15.4 2.4 GHz PHY layer. A review of low-rate wireless 
personal area network (LR-WPAN) solutions, including IEEE 802.15.4, can be found in [23]. 

The impact of interference in wireless sensor networks plays a very important role and can severely 
degrade the overall performance of the network and the efficiency of the upper layers. In our opinion, this 
aspect has not been sufficiently addressed in past years. In a dense sensor network deployment, where 
many nodes are periodically sending data to the sink, concurrent transmissions are highly probable. 
However, the probability of having a collision because of more than two concurrent transmissions 
is relatively low [18,20], thanks to CSMA-CA. In such conditions, the performance of the receiver 
depends on the overall amount of interferer signal energy and does not change with the number 
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of interferers [18,21]. Hence, solely the impact of one interferer on the capture probability will be 
considered in the mathematical analysis. 

On the other hand, we believe this study can also contribute to the identification of signal reception 
models for network simulators. In particular, as shown in Table 1, most common network simulators 
use signal-to-noise ratio threshold (SNRT)- or bit error rate (BER)-based models in order to decide the 
correct reception of a packet. 

Table 1. Signal reception models in network simulators [24]. SNRT, signal-to-noise ratio 
threshold; BER, bit error rate. 



Simulator 


GloMoSim 


ns-2 OPNET 


Signal Reception 


SNRT-based, BER-based 


SNRT-based BER-based 



In SNRT-based models, the packet is correctly received if the signal-to-noise ratio (SNR) is larger 
than a given threshold, whereas, in BER-based models, the packet reception decision is based on the 
BER, which is determined probabilistically depending on the value of the SNR. These models are 
rather simple, but have some drawbacks. In particular, SNR-based models do not take into account the 
impact of interference. This latter effect can be considered, in principle, by BER-based models, but the 
impact of the waveform of the interferer signals should be carefully considered, as it plays a significant 
role. Typically, conventional interference models are based on the assumption that the disturbance 
can be modeled as a Gaussian random variable; unfortunately, this is not the case of IEEE 802.15.4 
systems, where only a limited number of strong interferers is present. To counteract this problem, we 
mathematically analyze the impact of the waveform of the interferer on packet reception and obtain 
curves that are organized as specific look-up tables. Figures, such as those derived in Figures 4 and 6, 
can be used to provide accurate PHY models for network simulators. In that case, the conventional on-off 
behavior of SNRT-based models can be replaced by a probabilistic model, where the actual value of SIR 
leads to a given probability of packet loss. In other words, we provide a SIR-based signal reception 
model for the interference-dominant environments, where noise is not the serious cause of packet loss 
{i.e., enough transmit power is used or nodes are using the best links to reach the destinations in a dense 
sensor network deployment). Furthermore, Figure 7 shows that, in the case of non-coherent detection 
in an interferer-dominant environment, an on-off model can be also applied. In any case, behavior 
changes when thermal noise cannot be neglected. As a conclusion, the results of this paper on chip error 
rate (CER) and PRR (see Figures 6 and 7) can be used within network simulators in terms of look-up 
tables. That allows a fast characterization of the behavior of the PHY layer. 

The computational complexity of the model for the coherent detection is 0(1) (in big O notation). 
This makes it usable without intensive computational effort. For the non-coherent case, we show that the 
performance curve has a step-like behavior with the threshold at 0 dB. This simple model can capture 
the behavior of the non-coherent case without any computational effort. 

The rest of the paper is organized as follows: Section 2 describes CSMA-CA and the 2.4 GHz 
PHY of the IEEE 802.15.4 Standard. Sections 3 and 4 evaluate CER for coherent and non-coherent 
offset quadrature phase shift keying (0-QPSK), respectively. Section 5, obtains PRR in concurrent 
transmission by finding an upper bound, which transforms CER to data symbol error rate (DSER). 
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2. System Description 

Two different approaches are foreseen in IEEE 802. 15.4 to coordinate the data traffic: beacon-enabled 
mode and non-beacon-enabled mode [2]. In the former, periodic beacons are sent by the coordinator to 
synchronize the channel access, while in the latter, no synchronization is required. CSMA-CA is used 
in both modes, but in a slightly different way. In the beacon-enabled mode, every node synchronizes 
itself to the backoff slots determined by the coordinator; for this reason it is named slotted CSMA-CA. 
In non-beacon-enabled mode, every node has its own backoff slot timing, so that it is called unslotted 
CSMA-CA. Default values of the parameters and the number of sensing phases before assessing the 
channel idle are different in slotted and unslotted versions as a result of synchronization {i.e., two backoff 
slots in slotted, one backoff slot in unslotted). In both modes, each node waits for a random number of 
backoff slots, then channel sensing is performed. If the channel is found to be free, the node immediately 
starts the transmission. If the channel is found to be busy, the node turns back to the backoff state again. 
There is a maximum number of attempts a node can try to sense the channel. When this maximum is 
reached, the algorithm ends with a failure. Even though CSMA-CA avoids collisions by sensing the 
channel before transmitting a packet, a hidden terminal effect or sensing the channel as idle exactly at 
the same time by more than one node can cause a collision. As discussed in the previous Section, a 
collision does not always lead to a total loss. One of the packets, most probably the one with the better 
signal strength, can be successfully received. 

IEEE 802.15.4 PHY uses direct sequence spread spectrum (DSSS), and it operates over different ISM 
bands: the 868 MHz (EU) ISM band, the 915 MHz (US) ISM band and the 2.4 GHz global ISM band. 
Because of the global availability of the higher frequency band, currently transceiver manufacturers 
mostly have products working in this band {i.e., Chipcon CC2420, Freescale MCI 3224, etc.). In the 
paper, we focus on the 2.4-GHz frequency band, which utilizes half-sine shaped 0-QPSK modulation. 
In the PHY layer of the 2.4 GHz band, the signal is first modulated by forming a data symbol of four 
bits, and then, this symbol is mapped to 32-chip long sequences, obtaining 16 data symbols. The data 
symbol set modulated at the carrier frequency, /c, and carrier phase, (pc, can be written as: 



For each chip sequence, even-indexed chips are modulated to the in-phase (I) carrier, and odd-indexed 
chips are modulated to the quadrature-phase (Q) carrier of 0-QPSK. The offset between the I-phase and 
Q-phase is formed by delaying the Q-phase symbols one chip duration with respect to the I-phase. I{t) 
and Q{t) for each data symbol can be written as in [25]: 



sc.(t) = Ii{t) cos(27r/,t + 0c) - Qi{t) sin(27r/,t + 0c) ^ = 0, 1, 15 



(1) 



15 




(2) 



n=0 



15 




(3) 



n=0 
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where is ±1 with respect to the chip value and i is the index of the data symbols. The function, h{t), 
is half-sine shaped and can be expressed as: 

fcos(^) if -Tc<t<T^ 
I 0 otherwise 
Finally, the complex baseband signal is given as: 

scAt) = m + jQiity, 

the corresponding data symbol set in time domain will be expressed by using the complex envelope as: 

sciit) = ^[Scm exp (j(2vr/ct + 0c))] (5) 
where 3?[.] is the real part of the complex expression. 

3. Probability of Chip Error in Coherent O-QPSK 

In this section, we consider coherent communication and evaluate the chip error probability when 
there is one interferer with an unknown carrier phase and symbol timing. As the experiments 
in [16,18-20] suggest, we assume that the signal of the interferer is additive, and the probability of error 
is independent of the number of interferers. As a consequence, we will consider one interfering signal. 
Under such an assumption, the received O-QPSK modulated signal, in the time interval — Tc < t < Tc, 
can be written as: 



r(t) = ^[sc{t) exp (j(27r/,t + 0c))] +^[Si{t) exp (j(27r/,t + 



T,<t<T, (6) 



where sc(t) and sj{t) are the complex envelopes of the useful and interferer signals and (pc and 0/ are the 
carrier phases of the useful and interferer signals, respectively. There is no phase difference between the 
useful transmitter and the receiver, since we are assuming coherent demodulation. Therefore (f)c = 0, but 
the carrier phase of the interferer, 0j, is unknown. After some algebra. Equation (6) can be expressed as: 

r(t) = ^ [f (t) exp(j27r/,t)] - T, < t < T, (7) 

where the complex envelope is: 

r{t) = sc{t) + si{t) exp(j0,) -T,<t<T, (8) 

In the case of coherent demodulation, the I and Q components should be sampled at the instance 
corresponding to the maximum signal amplitude; hence for a sequence of chips, the I-phase should be 
sampled at the instants t/„ = 2nTc (see Equation (2)) and the Q-phase at the instants tg^ = {2n + l)Tc 
(see Equation (3)), where n = 0, 1, 2, . . .. This results in sampling the I and Q components of the useful 
signal at the indicated sampling instances in Figure la and Figure 3a. Due to the delayed sampling 
between the I and Q components, we will consider a Q-I constellation plane, where the Q component 
is delayed by Tc seconds with respect to the I component. We named the resulting constellation plane 
Q^-I to emphasize the delay. In the next two subsections, we derive the CER for the received signal in 
Equation (8), as a function of SIR, considering two cases: without and with pulse shaping. Thermal 
noise is neglected in the analysis, as we want to investigate solely the effect of interference on CER. The 
effect of noise will be considered in Section 5 to compare the analysis with the experimental results. 
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3.1. Without Pulse Shaping 



When we take the case of no pulse shaping into account, the shaping function, h{t), is a non-return to 
zero (NRZ) rectangular pulse, whose expression is: 



h{t) 



1 -T,<t<T, 
0 otherwise. 



(9) 



Accordingly, the complex envelope of the received signal in Equation (8) on the Q^-I constellation 
plane can be written as: 




^ ± ja/^ 1 h{t) exp{j<f)j) -T,<t<T,, 0 < 0/ < 27r (10) 



where C and / represent the energies of the useful and interferer signals, respectively. 0/ is the carrier 
phase of the interferer. Since the demodulation is coherent, perfect knowledge of the carrier phase 
and symbol timing is assumed for the useful signal. These parameters are supposed to be uniformly 
distributed random variables for the interfering signal. The pulse shaping function is constant throughout 
a chip duration; thus, random symbol timing of the interferer, t, does not give any contribution to the 
performance. On the contrary, the presence of an interferer carrier phase has a significant impact on the 
received signal, because (pi rotates the interfering signal on the Q^-I domain, as indicated in the second 
addend of Equation (10). The geometric representation of Equation (10) on the Q^-I constellation plane 
is depicted in Figure 1 . 

Figure 1. The impact of the interferer. (a) Useful signal; (b) Qd-^ domain representation 
of useful signal; (c) Interferer signal; (d) Q^-I domain representation of interferer signal; 
(e) Qd-l domain representation of received signal. 
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(c) 



(d) 



<|)i: Random carrier phase of interferer 



With respect to the transmitted symbol, the first addend in Equation (10) will be one of the 
constellation points on the Q^-I domain, as shown in Figure lb, while the second addend in Equation (10) 
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represents the interfering signal, which can be any of the points shown in Figure Id, but rotated according 
to the interferer carrier phase, 0/. The final sum in the equation can be interpreted as in Figure le; the 
sampling point of the received signal is located in one of the quadrants in the Q^-I plane, in agreement 
with the transmitted useful symbol, but, due to the presence of the interferer, it is shifted to a random 
location over the circle with radius \/l and the center at ^±y^, ±jy^j . The amplitude of the useful 
signal on the Q^-I plane is \/C, and the amplitude of the interferer signal is In the following 
parts, Pc and Pg will denote the probability of chip error and the probability of 0-QPSK symbol error, 
respectively. If we consider an ideal maximum likelihood demodulator [26], three distinct cases need to 
be investigated. 

Figure 2. Coherent Cases, (a) Vl < ^Jcj2\ (b) s/Cj2 < y/l < VC; (c) VC <y/l. 




(a) (b) (c) 



3.1.1. Case (a): < 

When vT is smaller than ^/Cj2, as shown in Figure 2a, the received signal never crosses other 
symbol regions. As a consequence, the demodulator will never make an error, and therefore: 



Pc = 0 
Ps = 0 

3.1.2. Case (b): < Vl < VC 

When vT is larger than (or equal to) a/C/2, the demodulator will start deciding erroneously, because 
of the stronger signal energy of the interfering component, which forces the received signal to be sampled 
in one of the adjacent quadrants on the constellation plane, as shown in Figure 2b. In such a condition, 
the probability of symbol error will be given by: 

p = ^ = 2— = — 

TT 27r TT 
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where 5*0 is the arc of the circle shown in Figure 2b, the angle, a, in radians, is equal to Sq/2, and a is: 



a = arctan 4 / ^ (11) 

V 

The sampling point of the received signal can only pass to the adjacent quadrants, and the receiver 
makes a one chip error per each transmitted symbol; as a result, will be the half of Pg'. 



^ An n 



3.1.3. Case(c): VC < Vl 



For increasing values of the circle will move to the other side of the origin (see Figure 2c). a can 
still be written as in Equation (11), but now: 



50 = a 

51 = 7r/2 — a 

52 = 2a - 7r/2 

When the received signal is in 5*2, the receiver makes two chip errors per each transmitted symbol. On 
the other hand, when it is in Sq or in ^i, the receiver makes one chip error per each transmitted symbol. 
Regarding the number of chip errors in the different quadrants, Pg and Pc can be written as: 

^ ^ 2(^0 + Si) + S2 1 
27r TT 4 

^ ^ 2(^0 + Si) l S2 

27T 2 27r vr 

The computational complexity of Ps and Pc is 0(1) in big O notation, since C and / are 
constant values. 

3.2. Half- Sine Pulse Shaping 

Using the pulse shaping given by Equation (4), the received signal in Equation (8) can be written, on 
the Qrf-I plane, as: 




-T, < t < Tc, 0 < 0/ < 27r (12) 



where the first addend represents the useful symbol and the second addend represents the interferer 
symbol. In the variable, t, in Equation (12), is the sampling instance for the interferer symbol, and 0/ is 
the carrier phase of the interferer, which shifts the pulse shaping function on the Q^-I domain. 

When there is perfect knowledge of the useful carrier phase and symbol timing, there is no difference, 
from the demodulator point of view, between half-sine-shaped or not shaped useful symbols, since 
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the sampling is done at the maximum amplitude instance. As a consequence, the first addends in 
Equation (10) and Equation (12) are equal; this can also be observed by comparing Figure lb with 
Figure 3b. When the pulse shaping function, cos(^ + 0/), of the interferer is included to the scenario, 
the interferer signal will be positioned on one of the diagonal lines shown in Figure 3d. As a consequence, 
the received signal plus interferer will be positioned over the dashed lines of Figure 3e. To obtain the 
error probability, we need to derive the probability distribution of the amplitude of the interferer at the 
sampling instance. 

Figure 3. The impact of of the interferer: the asynchronous case, (a) Useful signal; (b) Q^-I 
domain representation of useful signal; (c) Interferer signal; (d) Q^-I domain representation 
of interferer signal; (e) Q^-I domain representation of received signal. 




(c) (d) 



The random symbol sampling instance, t, and carrier phase, 0/, of the interferer can be modeled 
as uniformly distributed random variables. To simplify the notation in a pulse shaping function, we 
substitute the random variables as X = = cpj and Z = X + Y and denote the Probability Density 

Functions (PDFs) of X and Y as fx{x) and /y (y); then the PDF of Z, which is our interest, will be the 
convolution integral [27] given as: 



fz{z) 



fx{x)fY{z - x)dx 



' TT + 2Z 

47r2 


2 — ^ ^ 2 


1 

J 2tt 


n ^ ^ 3n 
2 — ^ ^ 2 


5tt-2z 
47r2 ' 


Stt ^ ^ Stt 
2 — ^ ^ 2 


0 
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Otherwise 


1 

2^ 


0< z <2iT 



(13) 



where: 
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I 0 otherwise I 0 otherwise 

The solution of the convolution integral in Equation (13) indicates that Z is a uniformly-distributed 
random variable in the interval 0 < 2; < 2n. Thus, the sampled point of the interferer signal is 
(i^/l ± jy'^) cos(z), and it always stays over the dashed lines in Figure 3d. If the interferer is 
transmitting ms or mo, the slope of the diagonal line is +1; on the other hand, if the interferer is 
transmitting mi or m2 then the slope is —1. The probability distribution of the positions of the interferer 
sampling point on the dashed lines are symmetric on the opposite sides of the origin, because of the 
cosine function. The amplitude in a quadrant is: 

i = ^/lcos{z) -'k/2< z <n/2 (14) 
Using Equation (14), the PDF of the amplitude is written as: 

fi{{} = , ^ 0<z<v^. (15) 



tta/I - U/VI 



Then, we can also obtain the probability that the interferer amplitude lies in an interval (a, b] using 
Equation (15): 

rv° 2 b 

Prob(a <t<b) = F,{b) - F,{a) = / ^ ^dt = -{sm-\ — ) - sm-\-^)) (16) 

where Fi{i) is the Cumulative Distribution Function (CDF) of /. 

Without loss of generality, as we did in Section 3.1, it can be assumed that the transmitter is 
transmitting 1113. Therefore, to obtain Pc and Pg, we should consider the following two cases: (a) 
Vl <VC, and (b) Vl>VC. 

3.2.1. Case (a): y/l < y/C 

When \/I is smaller than \/C, the demodulator never makes an error, since the received signal never 
crosses the other quadrants. Thus, Pc and P, are equal to zero: 



Ps = 0 



3.2.2. Case (b): y/l > VC 

When i, which refers to the amplitude of the interferer signal at the sampling instance, is larger than 
or equal to a/C/2, the receiver can decide for the wrong symbol with a probability of |. Regarding Pc, sl 
transition towards 11x2 or mi will result in one chip error per symbol, while a transition towards niQ will 
result in two chip errors per symbol. Therefore, the error rates when i > a/C/2 are: 
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111111 

P = + + - = - 

4 2 4 2 4 2 

In general, Pg and Pc can be written as a function of i as follows: 



0, i<VC 



i < Vc 

Vc <i<Vl 



Finally, we can derive the error probabilities when \/l is greater than \/C by evaluating the expected 



values of Ps{i) and Pc{i)'- 



p. = EiPM = Imd, = I ■ F,{Vc <'<^) = ll{l- s"r'(/f )) = ? - ^ si„-'(yf ) 

ft = EIPS)] = = I ■ F,(VC <i<^) = ll('^- Vj)) = ^ - ^ sm-'(\/f ) 

where /^(i) is given in Equation (15) and Fi(i) is given in Equation (16). 

The computational complexity of P^ and are 0(1) in big O notation, since C and I are 
constant values. 



3.3. Validation of the Analytical Model through Simulations 

The validity of our analytical framework has been tested using Monte-Carlo simulations. The 
simulation tool has been implemented in MATLAB (The Math Works Inc., Natick, MA, USA) using the 
complex envelopes of the signals. Asynchronous symbols are obtained by shifting the complex envelope 
in the time axis, while the random carrier phase is obtained through rotating the complex envelope on the 
complex plane. SIR values with a step of 0.3 dB from —10 dB to 5 dB are simulated. At each specific SIR 
value, 10, 000 random 0-QPSK symbols are generated for both the interferer and useful transmitter. The 
complex baseband representation of the 0-QPSK symbol is generated by using 100 points; therefore, 
the resolution of the time for the asynchronous 0-QPSK symbols was 10 ns. On the other hand, the 
resolution for the carrier phase was 27r/100. Random asynchronous interfering symbols with a random 
carrier phase are generated with uniformly-distributed parameters. The excellent agreement between 
simulation points and analytical curves is shown in Figure 4. Since the analysis in Sections 3.1 and 
3.2 is exact (no approximation is used), we expect a perfect overlapping between simulations and the 
analytical model. Asymptotically, there is a 3-dB difference between the results of half-sine shaped and 
non-shaped symbols. In half-sine pulse shaping, Pc goes to zero at 0 dB, while the threshold is 3 dB in 
the absence of pulse shaping. 
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Figure 4. Chip error rate: the case of coherent 0-QPSK demodulation. CER, chip error rate; 
SIR, signal-to-interference ratio. 
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4. Probability of Chip Error in Non-Coherent O-QPSK 

Normally, O-QPSK requires coherent detection [26]; however with the half-sine pulse shaping in 
2.4 GHz PHY of IEEE 802.15.4, the information-bearing part of the signal is not only the carrier phase, 
but also the complex-envelope of the signal. Thus, without recovering the carrier phase, just observing 
the phase changes of the complex envelope, transmitted information can be extracted. In 2.4 GHz PHY 
of IEEE 802.15.4, the phase of the complex envelope increases or decreases by an amount of 7r/2 every 
Tc seconds. In this sense, the behavior of the modulation is similar to that of minimum shift keying 
(MSK). The main difference between the two modulations is that they use different symbol mappings. 
In particular, symbol mapping in IEEE 802.15.4 is shown in Figure 5. 

Figure 5. Phase transitions of the complex envelope, (a) Example phase transition sequence; 
(b) All possible phase transitions. 
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Chip intervals in Figure 5a are marked as Tci,Tc2, ■■■■ In each even chip interval, Tc2jTcij a 
symbol (mo, mi, m2, ma) is transmitted, and in the odd chip intervals, Tc3,Tc5,..., there are phase 
changes (^4,^5,^6,^7)^ which can be considered as the transitions between two sequential 0-QPSK 
symbols. For instance, in Figure 5a, a sequence of phase changes are shown. At the beginning of 
the Tc2 chip interval {i.e., cut = 0), the phase of the signal is equal to the carrier phase, (pc', then at the 
end of the same chip interval, the phase reaches 0^ + 7r/2, which is indicated in Figure 5b by the arc, . 
In the next chip interval, Tcs, the phase reaches the 0c + value, which is represented by the arc, ts. The 
sequential chip intervals send the phase of the complex envelope, n/2, further or back with respect to the 
transmitted symbols. All possible phase transitions are shown in Figure 5b. This relationship between 
phase transition and transmitted symbol suggests the idea of constructing a non-coherent detector, which 
identifies the transmitted chip based on the phase transitions of the complex envelope. Using such 
a detector, will refer to the probability of the erroneous detection of phase transitions during a chip 
interval. In the absence of interference, the complex envelope of the received chip through a chip duration 
in the Q-I plane can be expressed as: 

R{ujt) = VCe^'-^'^'+'f'^\ wte[0,7r/2], 0ce[O,27r] (17) 

where C is the energy of the signal and (pc is the uniformly distributed carrier phase. The definition 
domain of cot is from zero to 7r/2, because the signal can only move 7r/2 further on the Q-I plane during 
a chip interval. The instantaneous phase of R{ut) in Equation (17) is: 

QRicot) = ±ujt + (j)c (18) 
and the phase of the received signal at the beginning and at the end of the chip interval is: 



Qr{0) = ct>c 
QRin/2) = ±7r/2 + 0^- 

Since the carrier phase can be eliminated by subtracting 6/j(7r/2) and Qr{0), the following discriminator 
demodulator can be applied [28]: 

en{7i/2) % Qr{0) (19) 

-j 

where the decision is based on the comparison of the phases of the complex envelope of the received 
signal at the beginning and at the end of the chip interval. The demodulator decides for a positive phase 
change, +j, if the phase at the end of the chip interval is higher than the beginning; otherwise, it decides 
for a negative phase change, — j. Firstly, we will analyze the case where the symbols of the interferer are 
synchronous (according to the chip interval) with the useful symbols. This assumption is not realistic, 
but provides the mathematical basis for the framework that will be used for the more realistic case, where 
the interferer and useful symbols are asynchronous. 
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4.1. Interferer and Useful Symbols are Synchronous 

When useful and interfering signals are synchronous, the received signal can be written as: 

R{ujt) = y^e^i^cMc) + y7eJ(fei-^t+</'/) (20) 

where tut = [0, vr/2], (pc = [0, 27r] and 0/ = [0, 27r]. be and bj represent the turning directions of the 
phase of the complex envelope of the useful and the interferer chips on the complex domain in Figure 5b: 
be, bi = +1 means a clockwise direction, whereas be, bj = —1 means the counter-clockwise direction. 
We denote Peg as the probability of erroneous detection when the phases of useful and interferer signals 
are turning to the same direction (be = &/), while Pdif denotes the probability of erroneous detection 
when phase transitions are toward the different directions (be 7^ bj). There are two cases; accordingly, 
Pcis: 

Pc = \Pe, + \Pd:f (21) 

Note that if the phase of the interferer signal turns toward the direction of the useful signal, the 
received signal in the Q-I plane will turn toward the direction of the useful signal, as well. As a 
consequence, P^q will always be zero. To obtain Pdif, we can consider the case where be = +1 and 
bi = —1; in fact, this is mathematically identical to the case of 6^ = — 1 and bj = +1. Motivated by this 
consideration, we can define the following success/error rule for the demodulator, based on the phase 
values at the beginning and at the end of the chip interval: 

6^(71/2) i Qr{0) (22) 

e 

where s means success and e means error. When be = +1 and 6/ = —1, the received signal is: 

R{ujt) = VCe^'(^*+'^^) + v^e^'(-'^*+'^^) (23) 

To obtain an analytical expression to ©^(wt), we assume that during a chip interval, R{ut) always stays 
on the definition domain of arctan (i.e., [—n/2 n/2]). Starting from Equation (23), the instantaneous 
phase function can be written as: 

[ ^/Csm{ut + (f)e) + Vl sm{-ut + 0,) ^ 

(c>R[ut) = arctan —= ^ (24) 

\ V C* cos{ujt + (pe) + yl cos(— cut + 0/) j 

With no loss of generality, we can assume (pc = 0. Now, after some trigonometric simplifications, the 
phases of the received signal at the beginning and at the end of the chip interval can be written as: 



0/?(O) = arctan 
0/j(7r/2) = arctan 



smicpj 



VC + Vl COS(0/) 

v^/)cos(0,) 



V/sin(0/) 

After replacing 9ij(0) and Qr{7i/2) in Equation (22), the decision rule becomes: 
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arctan I ^^cos(0J) \ g^j^^^g^j^ i ^' {4>i) i /25) 

' \/7sin(0j) )e \v^+y7cos(0/^ ' 

by recalling that arctan(-) is monotonically increasing in [— vr/l 7r/2], Equation (25) becomes: 



sm 0r 



(26) 



v^sin(0/) e cos 0r 

to keep the trajectory of R(ut) in the definition domain of arctan, we can restrict the definition domain 
of (J) J as [0, n/2]. With this assumption, \/C + vTcos(07) and y/l sm{(f)i) will be always positive. After 
some algebra. Equation (26) becomes: 

C-Icos\(j)j) I Ism\(f)j) 

e 

which gives: 

C I / (27) 

e 

It can be showed that (27) holds for the other possible domain definitions of R(cot); therefore, Pdif is: 



Pd^f = l C<I (28) 

Pd^f = 0 C>I (29) 



By replacing Equation (28) in Equation (21), the chip error probability becomes: 



Pc = l/2 C <I 

P, = 0 C>I 

4.2. Useful and Interfering Symbols are Asynchronous 

The condition in which interfering symbols can be asynchronous with respect to useful symbols 
increases the number of possible cases to be considered for the evaluation of QR^cut). In particular, 
the interferer signal can change its turning direction at a random instance, U, and the corresponding 
random angle is 0^ = utf. Consequently, the received signal during a chip interval can be written as: 

R{ut) = < (30) 

I ^/Ce^^^^^^+'f'^'^ + ^ti^^H-'^^'^) out = (0t , TT /2] 

where be defines the turning direction of the useful chip, while 6/^ and bi^ are the turning directions of 
the sequential chips of the interferer. In Equation (30), at the angle (pt, the interferer starts transmitting 
another chip. We denote this asynchronous interferer case probability as Pmix- The occurrences of Peg, 
Pdif (which are given in Section 4.1) and Pmix for the different combinations of be, bj^ and bj^ are shown 
in Table 2. 
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Table 2. Useful and Interferer Chip Phase Combinations. 





br 


br 


d 


Prohabilitv of Rrror 


+1 


+ 1 


+ 1 




Peg = 0 


+1 


+ 1 


— 1 


(Pl + 2(pt 




+1 


-1 


+ 1 




mix 


+1 


-1 


-1 


h 


Pdif 


-1 


+ 1 


+ 1 


h 


Pdif 


-1 


+ 1 


-1 




mix 


-1 


-1 


+ 1 


h - i<\>t 


P . 

^ mix 


-1 


-1 


-1 


h 


Pe, = 0 



By considering these occurrences, the chip error probability can be evaluated as: 

Pc = l/4Pe, + l/^Pdif + l/2Pmix (31) 

To obtain Pmix, we consider the case of 6c = +1.6/1 = +1» = d = (pj + 20^, which is 

mathematically identical to the all other Pmix cases in Table 2. Following the same approach followed 
in Section 4.1, we can write the decision rule as: 

VC-Vl cos(0/ + 2(j)t) I Vl sin(0/) 

^= < — ^= ^= i-'^J 

VI Sin(07 + 2(f)t) e y/C +VI COS(07) 

which can be simplified to: 

C + VCI{cos{(pi) - cos(0/ + 2(pt)) - / cos(20t) ''^ 0 (33) 

e/s 

where the sign in the inequality depends on the sign of denominators in Equation (32). Equation (33) 
can be rearranged as: 

C/I + >/cJl{cos{(f)i) - cos(0/ + 2(pt)) - cos(20t) 0 (34) 

e/s 

Pmix is obtained by using the test in Equation (34) considering the different signs of the denominators 
in Equation (32) and the definition domain of the arctan function. 

4.3. Validation of the Analytical Model through Simulations 

Simulation results for the cases of synchronous and asynchronous useful and interfering symbols are 
shown in Figure 6. Since the analysis in Section 4 is exact, we expect a perfect agreement between 
simulations and analytical model. The results of Figure 6 confirm this intuition. 



Sensors 2014, 14 



5638 



Figure 6. Non-coherent chip error rate. 
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5. Probability of Data Error in IEEE 802.15.4 

The received path of the commercial 802.15.4 compatible transceivers, such as TICC2420 [29] and 
Freescale MC13224 [30], is composed of two main blocks: an 0-QPSK demodulator and a de-spreader 
to obtain the data symbols from the received chips. In the previous section, we have obtained the chip 
error rate, Pc, by assuming both coherent and non-coherent detection; this error probability is used here 
to derive bounds for the symbol error probability, Pd, and for the corresponding packet reception rate, Pr. 

5. 1. Chip Error Rate to Data Symbol Error Rate 

As discussed previously, in the 2.4 GHz PHY of IEEE 802.15.4, data symbols are mapped to the chips 
after the demodulation. Because of the symmetry in the data symbol set used by the 802.15.4 2.4 GHz 
PHY layer, the error probability of the data symbol, Sq, is equal to the symbol error probability, P^ [31]: 

Pd = Pie\So) = P{e\S,) for 2 = 1, . . . , 15 (35) 

Hence, with no loss of generality to obtain an analytical expression for Pd, we can find the union upper 
bound of P(e|5'o) expressed as: 

15 

P{e\So)<J2PiS^\So) (36) 
1=1 

where P(Si\So) represents the probability of deciding for the symbol. Si, once symbol 5*0 has been 
transmitted. Each probability in the summation of Equation (36) can be formulated as a function of 
Hamming distance, hi, between 5*0 and Si, and the chip error rate, Pc, as: 



P{Si\So) = P{h,,P, 



k&P^'^'i^ - Pc)'''^' + ^t^^ (!^)Pc{l-Pc)''-' if /i. is even 



E Ln^ (!^)P^il-Pc)''-' if /l, is Odd 

*~ 2 

(37) 

By substituting Equation (37) in Equation (36), we obtain an upper bound that can be applied to any 
symmetric symbol set. The bound will be used in the next sub-sections to obtain expressions for the data 
symbol error probability of both coherent and non-coherent detection. 
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5.1.1. Coherent Detection 

In the case of coherent detection, we can calculate the Hamming distance of each chip sequence. In 
particular, the recurrences of the Hamming distances of the data symbols is given in Table 3. 

Table 3. Recurrences of unique Hamming distance values. 



Hi 12 14 16 18 20 

Reoccurrence (i?/) 2 2 3 2 6 



Using these recurrences, the symbol error probability can be bounded as follows: 

N 

Pd<Yl = Pe) + 2P(14, P,) + 3P(16, PJ + 2P(18, PJ + 6P(20, PJ (38) 

1=1 

where P{hi, Pc) is taken from Equation (37) and hi and are given in Table 3. z is the index for columns, 
and = 5 is the number of columns in the table. 

Note that, for a fixed value of Pc, the computational complexity of the previous upper bound is 0(1). 

5.1.2. Non-Coherent Detection 

In the non-coherent case, we can only identify the phase transitions during the chip intervals. For a 
32-chip long sequence, there will be 31 phase transitions, as shown in Table 4. 

Table 4. Phase transitions of chip sequences in 2450-MHz PHY of IEEE 802.15.4. 



Data Symbol 


Data Bits 


Phase Transitions 


5*0 


0000 


+ + + + + . + + + + . + . + + + .. + + . + + .. 




1000 


+ .. + + + + + + . + + + + . + . + + + .. + + . 


^2 


0100 




^3 


1100 


+ + -- + + - + + -- + + + + + + . + + + + . + . 


^4 


0010 


. + . + + + .. + + . + + .. + + + + + + . + + + 




1010 





s. 


0110 






1110 






0001 


.. + + + + + + — + + . + — + + .. + .. + + 


S9 


1001 


— — -1- — 




0101 


— — -1- — -1- 


Su 


1101 


— — -1- 


S12 


0011 


+ . + — + + .. + .. + + — + + + + + + — + — 




1011 


+ - + — + + .. + .. + + — + + + + + + — 


Si4 


0111 


— + + . + — + + .. + .. + + — + + + + + 


Sl5 


nil 


+ + + + — + + . + — + + .. + .. + + — + 
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Where + means a 7r/2 increase in the phase and — means a 7r/2 decrease in the phase. We can think of 
the phase transition sequence of a data symbol in Table 4 as a vector in Hamming space and calculate the 
Hamming distances between these vectors. The recurrences of Hamming distances from a data symbol 
to the other data symbols in the set are given in Table 5. 

Table 5. Recurrences of unique Hamming distance values. 



hi 13 14 15 16 17 18 31 

reoccurrence (r^) 2 2 3 3 2 2 1 



The union upper bound is obtained as: 

N 

Pd<Y, ^*^(^'' = 2^(13, Pc) + 2P(14, P,) + 3P(15, P,) + 3P(16, P,) + 2P(17, P,) 

i=l 

+ 2P(18,Pe) + P(31,P,) (39) 

where hi and are given in Table 5. i is the index for columns, and = 7 is the number of columns in 
the table. 

Again, the computational complexity of the upper bound is 0(1). 

5.2. Packet Reception Rate 

In 802.15.4 2.4 GHz PHY, four data bits are mapped to a data symbol [2]. When a packet of 6-bytes 
is transmitted, the probability of successfully receiving the packet, P^, is calculated as: 

Pr = {l- PdT (40) 

It is worth noting that the upper bounds on in Equation (38) and Equation (39) will become lower 
bounds in Equation (40). Pr as a function of SIR is shown in Figure 7 for the coherent and non-coherent 
demodulation schemes considered in Sections 3 and 4. Comparison between analysis and simulations 
is also shown. Results for the case of no thermal noise (cross-points for coherent and filled circles for 
non-coherent) have been obtained for a PHY service data unit (PSDU) of 14 bytes. The figure shows 
that the agreement between analysis and simulations is very good. 

Experimental results shown in Figure 7 of [19] suggest 7 dB for the threshold of SNR, where the 
performance of the system, in terms of PRR, completely switches from zero to one. In this sense, the 
interval from SNR = 7 dB to SNR = oo should be adequate for representing a great variety of different 
SNR operation conditions, where the noise is not dominant, but the performance of the system is mostly 
determined by the interferer. This interval results in a considerably narrow region of PRR (the shaded 
area in Figure 7), and the experiment results in [16,18-20] fall into this region; therefore, it would not 
be wrong to conclude that the mathematical analysis presented here captures the essential features of 
the real- world performance, where the dominant factor of performance degradation is interference. Two 
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different values of PSDU are considered in Figure 7: 14 bytes (experiments in [20]) and 127 bytes (the 
maximum allowable, which is used in experiments in [16,18,19]). Finally, It is also worth noting that 
0 dB is the absolute minimum SIR necessary to successfully detect the useful packets in non-coherent 
demodulation. 

Figure 7. Packet reception rate (PRR). 
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6. Conclusion 



The effect of concurrent transmission on the performance, in terms of chip, symbol and packet error 
probability, of IEEE 802.15.4 systems was investigated. To this aim, an analytical model, which starts 
from the PHY characterization of the IEEE 802.15.4 standard, was proposed, and the results were 
validated using simulations. Since the system performance does not depend on the number of interfering 
devices, but rather, on the overall amount of interfering power [18,20], the case of one single interferer 
was considered. Both non-coherent and coherent demodulation were analyzed in the absence of thermal 
noise. The joint effect of thermal noise and interference was taken into account by means of simulations. 
For the case of non-coherent demodulation, we found that 0 dB is the absolute minimum necessary 
value of SIR for the successful detection of the useful packets. We have also shown that the proposed 
analysis is in agreement with the measurement results of the literature [16,1 8-20] under realistic working 
conditions. We conclude that the mathematical framework can be used to provide analytical explanations 
to the experimental results shown in [16,18-20]. 
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